The effect of confinement on the interaction between two like-charged rods 
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Monte Carlo simulations were employed to study two charged rods confined between two unlike 
charged plates, all immersed in a model electrolyte. Recently, it was shown that two rods immersed 
in a multivalent counterion solution may show attraction (PRL 78, 2477 (1997)). Here we show 
for a monovalent electrolyte that rod-rod attraction and repulsion can switch sign depending on 
confinement and ionic size. We also propose a simple self-assembling mechanism which may be 
helpful to understand the DNA-lipid bilayers complexation. 
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Confinement effects appear in a wide range of length 
scales. They are observed from very large scales up 
to subatomic scales Q, passing though molecular levels 
[H, El]. At molecular levels, confinement induces order- 
ing, produces density fluctuations, and provokes dielec- 
tric changes on liquids, and even enhances miscibility in 
polymer blends [3|-l5j. In colloidal science, it is known 
to produce attractions that have yet to be explained [6J, 
and in biological setups, the presence of plate-like, paral- 
lel, charged lipid bilayers was shown to induce attraction 
between like-charged DNA molecules (7- 

In soft-matter and at nanoscale levels, confinement af- 
fects many body electrostatic forces, short range correla- 
tions, and thus, macroparticles double layer interactions. 
Additionally, other associated phenomena such as coun- 
terion release appear, as a function of the distance from 
the macroparticles to the approaching, unlike charged, 
confining surfaces. Counterion release is understood as 
the expulsion of counterions from the space between the 
macroparticles and the approaching surfaces 04131 , and 
it is claimed to be the driving force for higher-order self 
assembling 0, EH- This conclusion was built from 
DNA lipid bilayers complexation experiments, assuming 
that the cationic lipid bilayer fully neutralize the DNA 
charge. As we will show, this condition is not fulfilled. 
For the lamellar phase, this system may be modeled by a 
couple of negatively charged rods (DNA) , sandwiched by 
parallel positive charged plates (bilayers). Nevertheless, 
aside from important contributions dealing with paral- 



lel rigid polyelectrolytes [13l - ll8l |. there are no theoretical 



works on the effects that confinement may produce on 
their interaction. Thus, we perform Monte Carlo (MC) 
simulations of this system for elucidating these effects 
and the associated phenomena. 

NVT MC simulations were performed for two in- 
finitely large, hard, parallel, negative rods confined be- 
tween two hard, parallel, positive plates of finite thick- 
ness, which are immersed in a 1-1 restricted primitive 
model (RPM) electrolyte, unless otherwise specified (see 
Fig. [1]). RPM consists of hard spheres with a centered 
point charge, such that their electrostatic interaction is 



U E (rij)-- 



'-, being /3=l/fcs? 1 , Ub the Boltzmann con- 

and Zj the 



stant, T=298K the absolute temperature, 
valences of sites i and j, £b — ^-—7-1AA the Bjerrum 
length, e=78.5 the dielectric constant, and the inter- 
particle distance. For simplicity, plates, rods, and solu- 
tion have the same dielectric constant. The simulation 
box has sides lengths L X =L Z = 200 A, and L y = 125 A. 
The origin of coordinates is set at the box center. Plates 
are located parallel and symmetric to the z=0-plane, 
with a surface-surface separation distance r. Rods axes 
are placed parallel to the y axis at x=±(R+t/2), z=0, be- 
ing R the rods radii and t the rods surface distance. We 
assigned to plates a charge density of co= 0.229C/m 2 
(+e per 70 A ) for each surface, in correspondence 
with a fully charged cationic lipidic membrane [19(. We 
place +e-charges over a triangular grid, which is ran- 
domly moved in the x and y directions to mimic a con- 
tinuous charge distribution. Rods have an — e-site each 
1.7 A, as DNA molecules [3]. Rods have a radius of 
i?=10.5 A, consistent with a hydrated DNA molecule. 
Plates are 5 A thick to allow fluid-fluid correlations be- 
tween the interplate and bulk sides of the plates [2fj| . 
These correlations are expected in real 40 A thick bi- 
layers, due to their low dielectric constant. We fixed 
the concentration of the 1-1 electrolyte to p s =0.1 M, 
and assigned its diameter to a=4.25 A, unless other- 
wise indicated. Additional anions are added to make 
the system electroneutral. Periodical boundary condi- 
tions are set for all directions. The Ewald summation 
formalism is employed to deal with Coulomb interactions 
2ll 22| . To access all phase space volume, i. e., to allow 



for ion interchange between the confined and unconfmcd 
regions, movements having a large maximum displace- 
ment are done. Electrostatic contributions to the forces 
acting on the macroparticles (rods and plates) are ob- 
tained by F e i={J2iJ2j -W E (rij)), where i runs over 
the sites of the reference macroparticle and j runs over 
all other sites. The contact contribution is obtained 
by integrating the ions contact density,p(s=cte), F c = 
— J ksT ' p(s=cte)nds , where s refers to the macropar- 
ticle surface and n is a unit normal vector. These two 
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FIG. 1. Scheme of the system showing regions ri, r2, and r-j,. 
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FIG. 3. Negative of p e i(x,z) for h/a = 2.5 (a), 1.5 (b), 1.0 
(c), and (d). Blue indicates a larger concentration of an- 
ions than cations and red the opposite. Darker tones indicate 
larger absolute values of p e i(x,z). The inset of each figure 
schematizes the charge distribution. 



FIG. 2. a) Cation, o, and anion, •, densities for region n. 
b) Normalized charge densities for regions n (•), ri (o), and 
(rs) x. c) Normalized charge density of region | z \ < r/2. 



contributions to the force are interdependent. 

First simulation experiment consists of studying the 
system as a function of h=(r — 2R)/2, for t=0, being 
h the rod-plate surfaces distance. To study the ion re- 
lease mechanism we focused on the following interplate 
(| z | < r/2) regions: (n) \x | < 2R, (r 2 ) 2R < | x \ 
< 2R+3l D , and (r 3 ) 2R+3l D < | x | (see Fig. [J), be- 
ing l D ={e/{2f3e 2 p s )} 1 / 2 =9.6A the Debye length. The 
average density of cations and anions, pi— J pi{v)dv, 
are presented in Fig. [5] a). Figure [5] b) shows the nor- 
malized charge density, o rj / (2<jq), for region rj, where 
cr rj =-fi- J r . Pei{r)dr, Aj is the area of the j region pro- 
jected on the plate, and p e /(r) is the volume charge 
density including anions, cations, and rods. Figure [2] 
c) shows the normalized charge density of the inter- 
layer region, a in /(2a ), being a in = ^X^ ^jOY, and 
At=J2jAj- In Fig. [2]b) is seen that in r 3 the induced 
charge practically compensates that on the plates for all 

h. This is something expected since this region is rel- 
atively far from the rods, and plates are relatively far 
from each other: r>2lA«5a. On the other hand, lo- 
cal electroneutrality is violated for the other two regions, 

i. e., the charge of the confined fluid does not compen- 
sate the corresponding charge of the plates [231 ] . It is 
seen that the charge of region r\ always overcompen- 
sates it, whereas region always undercompensates it. 
Both effects, in turn, induce charge polarization in the 
fluid outside the plates. However, the overall charge of 
region r\ + r 2 does not compensate the charge of plates, 



but produce violations of charge electroneutrality. This 
yields the <Ji n /{2ao) plot shown in Fig. H]c). This ef- 
fect was predicted by integral equation studies for con- 
fined electrolytes. Furthermore, it was also shown that a 
correlation exists between the inside and outside of the 
plates [2Gj . Here, we found both effects to be in agree- 
ment with those results (see also Fig. . There is under- 
compensation of the plate's charge for h<0.5a and over- 
compensation for h>0.5a (Fig. H]c)). Electroneutrality 
is yield for large h values. It is also seen from Fig. [2] 
a) that ion release starts at h « 4a, as plates approach 
each other. However, cation release ends at h=1.5a since 
region r\ runs out of them. On the contrary, anion re- 
lease is strongly produced for h<1.5a. Hence, both facts, 
running out of cations and strongly releasing anions for 
/i<1.5a produce the decrease of negative charge at region 
n (solid symbols of Fig. [2]b)), which in turn manifest as 
the maximum of <Ji n /(2ao) in Fig.[2]c). Electroneutrality 
is violated at the expenses of keeping a constant chemical 
potential 2^, 23[. As we discuss below, Oi n determines 
the electrical contributions to the force on plates. 

The charge distribution profile around rods and plates 
is defined by p e i{x, z)=z + g + {x, z; r) + z_5_(x, z; t), be- 
ing g+(x, z; t) (g_(x, z; r)) the cation (anion) distribu- 
tion profile, given t; p e i{x,z) is shown in Fig. [3] for dif- 
ferent values of h. For large h, the plates double layers 
are uncorrelated from those of the rods and from those 
outside the plates (not shown). For smaller h, the dou- 
ble layers start to interact, becoming strongly correlated. 
This is already seen in Fig.[3]a), where the positive charge 
around the rods distorts due to the presence of the plates 
double layers and viceversa. Consequently, depiction re- 
gions appear between rods and plates. As pointed out 
above, from a mechanistic point of view, these depletion 
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FIG. 4. a) Forces on rods per unit of rod area (pressure) as a 
function of h for rod-rod contact. Crosses were obtained when 
no plates were included in the simulation box. b) Forces on 
plates per unit of plate area (pressure) as a function of h. c) 
Pressure on rods as a function of t for different h. 

regions are linked to partial release of cations and anions 
from the interplate region towards the bulk, thus, due to 
the plates and rods double layers interactions, i. e., the 
particles charge and size correlations. However, from a 
fundamental point of view, the depletion regions are due 
to the need of the system to increase particles accessi- 
ble volume and, hence, the entropy, and ultimately, due 
to the conservation of the chemical potential. Figure [3] 
b) is built for hml.5a, where cations were practically ex- 
pelled out from the interplate region. There, small red 
peaks recall us the presence of positive ions. Additionally, 
the anionic depletion regions were largely increased. For 
smaller h, no more cationic peaks are observed (Fig.[3]c) 
and Fig.[3]d)), and the anionic depletion regions enlarge. 
This is just a consequence of excluded volume, since no 
ions can enter the region in between rods and plates. In 
addition, the regions at x—0,z—±y/a(R + a/4), where 
cations used to place, are now populated by anions. The 
insets in Fig. [^schematize the ionic distribution for these 
configurations. 

Effective rod-rod forces per unit area (pressures) are 
shown in Fig. 2] Note that system symmetry makes the 
rod-rod force to have only components on x. This com- 
ponent is shown in Fig. S]a), and a positive sign means 
repulsion. Similarly, just the plate-plate z component of 
the force, per unit area, is shown in Fig. H]b). Figure 
S] a) shows that the influence of plates on rods force is 
non negligible only for h<4a. For larger separations, the 
obtained force is equal to that obtained when no plates 
are included in the system (crosses). This is a conse- 
quence of the fact that plates and rods double layers are 
uncorrelated, as explained above. This force is composed 
by a large electric repulsive contribution plus a positive 
contact contribution, which arises as a consequence of the 
large adsorption of counterions at x—0,z—± \J a(R + a/4) 



[18| . The contact contribution starts decreasing for /i«4a 
and reaches zero for /i«1.5a. This decreasing contact 
contribution is a consequence of the hight decrease of the 
cations peaks at x—0,z— a(R + a/4). The attractive 
contribution of this force for h<a is due to the anion 
peaks that grow at the outer sides of the rods, which 
overcompensate the contribution of the anion peaks at 
x=0,z= ±y/a(R + a/4). This is seen in Fig. [3] c) and 
[3] d) . Finally, the decrease of the electric contribution 
is due to the increase of the anionic concentration in T2 
and r3, which partially compensates the rod- rod non- 
screened electric force. Both contributions produce an 
ever decreasing total force, which yields zero for rods- 
plates contact, h=Q. 

Figure 2] b) shows negative values of the total plate- 
plate force for all distances, indicating that plates attract 
each other. For large values of h both contributions to 
the total force tend to vanish due to the independence of 
the double layers. As h is decreased the contributions be- 
come negative. The electric contribution, however, shows 
an absolute minimum at /i«1.5a. Thus, for h<1.5a it in- 
creases yielding a positive value for h=0. These features 
recall us the h dependence of ai n /{2a ) (see Fig. H]c)). 
In fact, both curves are practically equally shaped, in- 
dicating that the electrical component of plate pressure 
is given by P e i~(a 2 out - of„), where a out is the induced 
charge density for the unconfined fluid [2(| ■ Hence, coun- 
terion release, which is related to the charge density, in 
turn is linked to the electrical component of the pres- 
sure on plates. On the other hand, depletion regions are 
the cause for the attractive contact contribution to the 
plates forces. As was shown in Fig. 02 depletion regions 
enlarge as h decreases, explaining the general increase of 
the contact attraction. It should be noted that it is the 
contact force the main contribution to the net force on 
plates. In fact, for rod-plate contact, this is the com- 
ponent that yields attraction, in contrast to the electric 
contribution. Finally, Fig. @]c) shows the forces on rods 
per unit area as a function of the rods surface distance, t, 
for different degrees of confinement. It is clearly seen how 
confinement affects the effective force, producing practi- 
cally null values for all t and h—0. In fact, the work 
per charged site to approach the rods from t=4a to t=0 
is 2.13/cbT '/ site for unconfined conditions, Q.TSksT / site 
for h=a, and 0.02k B T/site for h=0. 

For the parameters employed, experiments show the 
spontaneous self assembling of bilayers (plates) and poly- 
electrolytes (rods) into a compact structure [10j| . Our 
results clearly indicate that plates squeeze the rods. On 
the other hand, although our model shows that the neces- 
sary energy to pack the rods is practically zero for highly 
confined conditions, it does not predict attraction. Nev- 
ertheless, there is a very important difference between 
real experiments and our model; that is the employed 
rod concentration. Experimentalists need a sufficiently 
large polyelectrolyte concentration to obtain the compact 
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FIG. 5. Force acting on rods as a function of the ionic size. 



structure: As more polyelectrolyte enters the midplane, 
the accessible volume fraction decreases, enhancing short 
range correlations, which in turn produce higher adsorp- 



tion at surfaces [24j. This effect may produce a net at- 



tractive forces on rods. In order to show this, we did two 
final computer experiments for h=t=0. A): we studied 
the influence of the ion size, a. B): we add a macroion 
species, with valence —5 and diameter 13.4 A, at the iso- 
electric point concentration (when the membranes charge 
equals the macroions charge) [Io| . Here, we also include 
all plates counterfoils (valence —1, diameter 4.25 A) and 
the necessary cations (valence 1, diameter 4.25 A) to keep 
the system electroneutral. We wish to mimic a DNA-lipid 
complex counting on a finite DNA concentration, i. e., 
we expect the macroions to behave as monomers of the 
DNA molecule. Hence, the macroions size-charge rela- 
tion was chosen in such a way that the electric field at 
their surface equals that at the rods surface. 

Results showing the influence of the ion size (A) are 
in Fig. [5] It is observed that the total force becomes at- 
tractive for a>4.25A. More importantly, it is clear the 
influence of the ionic size on the rod-rod effective inter- 
action. That is, the larger the little ions, the lower the 
accessible volume and so, the system entropy is forcedly 
reduced. Consequently the system reacts pushing more 
ions against the plates and rods surfaces, to gain acces- 
sible volume and hence, entropy. From a mechanistic 
point of view, this translates into a larger contact at- 
tractive force. This is seen in Fig. [5l where the con- 
tact contribution clearly rules the total force for a large 
enough ionic size. In the second experiment (B), we ob- 
serve (not shown) the macroions to replace the plates 
countcrions and to crowd inside the interplate region and 
around the rods. Thus, the rods induce macroions lines 
next to them, which, we argue, would play the role of a 
finite rods concentration. We obtained an effective rod- 
rod pressure of — 0.81MPa, and a plate-plate pressure of 
— 10.38MPa. Both are attractive pressures that signature 
the self-assembling tendency of these systems. Hence, 
these results suggest that complexation into highly pack- 
aged structures should be obtained for a certain concen- 



tration of large macroions (polyelectrolytes) . This is ex- 
actly what experiments show, where small ions leave the 
confined region in favor of larger ones, to increase the 
system entropy. 

In brief, not only handling valence and concentration 
of added salts ma y r esult in attractive forces between 
macroparticles ll5Ml8l ] but also confinement may pro- 
duce this effect. It was shown how confinement of rods by 
opositely charged plates reduces their inherent repulsive 
effective interaction forces as plates approach each other. 
This is a consequence of the need of the system to lower 
the free energy by increasing entropy. Nevertheless, the 
process is linked to several interesting phenomena such 
as counterion release, local violations of clcctroncutrality, 
and to the appearence of ionic depletion regions. Here we 
show that attraction between rods and plates depends on 
a delicate balance between charge and excluded volume 
of all particles in the system. 
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